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Abstract 

We consider neutron radiative /3-decay, n — t- pe'v'ej, aud compute the T-odd momentum corre- 
lation in the decay rate characterized by the kinematical variable ^ = • (le x k) arising from elec- 
tromagnetic final-state interactions in the Standard Model. Our expression for the corresponding 
T-odd asymmetry is exact in 0{a) up to terms of recoil order, and we evaluate it numerically 
under various kinematic conditions. Noting the universality of the V-A law in the absence of recoil 
order terms, we retain the parametric dependence on masses and coupling constants throughout, 
so that our results serve as a template for the computation of A^^^ in allowed nuclear radiative 
/3-decays and hyperon radiative /3-decays as well. We consider the pattern of A^^ which appears 
^ in nuclear decays and show that A^^ can be suppressed in particular cases, facilitating searches 

for new sources of CP-violation in such processes. 
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I. INTRODUCTION 



Radiative /3-decay offers the opportunity of studying T-odd momentum correlations wliicfi 
do not appear in ordinary /3-decay [1]. We consider a correlation cfiaracterized by tfie kine- 
matical variable C, = 1;^ • (le x k), so tliat it is botfi parity P- and naively time-reversal T-odd 
but independent of the particle spin. Its spin independence renders it distinct from searches 
for permanent electric-dipole-moments (EDMs) of neutrons and nuclei. The inability of 
the Standard Model (SM) to explain the cosmic baryon asymmetry prompts the search for 
sources of CP-violation which do not appear within it and which are not constrained by 
other experiments. A triple momentum correlation in radiative /3-decay is one such example, 
as we shall illustrate; under the CPT-theorem, T-violation is linked to CP-violation. A 
decay correlation, however, can be, by its very nature, only "naively" or "pseudo" T-odd, 
that is, only motion-reversal odd. As a result, although the appearance of a T-odd decay 
correlation can be engendered by sources of CP- violation beyond the Standard Model, it can 
also be generated without fundamental T- or CP-violation. In this paper we compute the 
size of the T-odd momentum correlation in radiative /3-decay simulated by electromagnetic 
final-state interactions in the SM [2J. This is crucial to establishing a baseline in the search 
for new sources of CP-violation in such processes. Our work is motivated in large part by 
the determination that pseudo- Cher n- Simons terms appear in SU(2)lxU(1) gauge theories 
at low energies - and that they can impact low-energy weak radiative processes involving 
baryons [3H3]. In the SM such pseudo-Chern-Simons interactions are CP-conserving, but 
considered broadly they are not, so that searching for the P- and T-odd effects that CP- 
violating interactions of pseudo-Chern-Simons form would engender offers a new window on 
physics beyond the SM [6]. 

Searches for T-violating decay correlations in neutron and nuclear /3-decay have a long 
history. The best experimental limits are on the so-called D term, which appears as the 
triple correlation DS ■ (le x 1^,), where S is the polarization of the decaying particle [71 [8]. 
These limits still greatly exceed the size of the D correlation expected from SM final-state 
interactions [HI HO]- Radiative /3-decay offers the possibility of forming a T-odd correlation 
from momenta alone; to our knowledge such a possibility was first considered in the context 
of Ki^^ decay [TT] . The T-odd asymmetry computed in Ref . [TT] from electromagnetic final- 
state interactions has recently been recalculated and is in significant disagreement with the 
earlier result [12j . 

In this paper we evaluate the T-odd asymmetry in radiative /3-decay from electromag- 
netic radiative corrections in the SM and focus on the neutron case: — )■ p{pp)+e~(le) + 
^e{^v) The motion-reversal-odd terms in the decay rate, which mimic the appearance 

of T-violation, are engendered by the interference of the tree-level amplitude with the imagi- 
nary part of the 0{a) corrected amplitude, which is determined by the physical two-particle 
cuts and hence mediated by the scattering of particles on their mass shells [HI HSl IE]- In 
what follows, we detail the computation of the interference terms and their components, as 
well as the resulting numerical integration over the allowed phase space to yield the T-odd 
asymmetry A^. Our results are exact in 0{a) up to corrections of recoil order, namely, up 
to terms of 0{e/M), where e is an energy scale which is small with respect to the nucleon 
mass M . This certitude is guaranteed by the small Q- value of the decay, so that e ^ M, 
and by Low's theorem [15]. The natural scale of hadron excitations is set by the pion mass 
m^r; consequently, in neutron radiative /3-decay e <^ 771^^ as well, and non-electromagnetic 
final-state interactions cannot contribute to the physical two-particle cuts. This is in con- 
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FIG. 1: Contributions to n{pn) — )• p{pp) + e~{le) + Ve{lv) + up to corrections of recoil order. 
The effective weak vertex is denoted by ® and is controlled by the Fermi constant Gp. The diagram 
enumeration is utilized in our calculation of the T-odd asymmetry. 



tradistinction to Kj^^ decay for which such contributions are appreciable, albeit relatively 
small [16]. We relegate intermediate results essential for our final results but yet nonessential 
to the flow of our discussion to Appendices. Since we neglect all terms of recoil order, our 
results can be immediately translated to the computation of in nuclear and hyperon 
radiative /3-decays as well. We consider, in particular, the pattern in predicated by 
variations in the Q-value of the decay as well as by the nuclear dependence of the axial- 
vector coupling constant. We assess the size of undetermined corrections before offering a 
final summary of our results. 



II. FORMALISM 



We work in a simultaneous expansion in the electromagnetic coupling constant e and in 
e/M, so that the leading contributions to neutron radiative /3-decay are from the diagrams in 
Fig. [Tj At this order the baryons are effectively structureless, and the contributions arise from 
bremsstrahlung off the charged particle legs of ordinary /3-decay, yielding a gauge-invariant 
result [E]- Employing the notation and conventions of Ref. [17], the decay amplitude is: 

A^tree = A^Ol + Mo2 , (l) 

with [H [E] 

MQl{le, k,pp) = ^^^^^ Ue{le) ^^'' ^ 1 - l^)Vv{lv)Up{j)p)Y " \l^)Un{Vn) , (2) 

A402(/e, k-.Vp) = - ^^^^ ^e(^e)7p(l " l5)Vv{lu)Up{pp) ^^^ ^ 7^(1 - \^^)Un{Vn) , (3) 

where is the photon polarization vector and A = gAlQvi noting that gy and qa are the 
vector and axial-vector weak coupling constants of the nucleon, respectively. We explicitly 
include the arguments in the momenta Ze, k, and Pp for later convenience. 

The branching ratio and photon energy spectrum for this process have been computed 
previously [18] [T9] . The expressions which follow from Eq. ^ are consistent with the 
experimental results [201 EI]- The next-to-leading order terms in the small-scale expansion. 
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i.e., those of 0{e/M), have been computed in heavy-baryon chiral perturbation theory and 
are no larger than 0(0.5%) of the leading-order result |19j - this is some twenty times smaller 
than the current experimental sensitivity [21j. In what follows we neglect all recoil-order 
terms and consider the 0{a) corrections to the amplitude of Eq. For future reference, 
employing lepton and hadron tensors, we note that 



where M„, Mp, and u refer to the neutron mass, the proton mass, and the photon energy, 
respectively, and 

^ee^ps ^ _Q4M^Mp {ml -h-k) ((1 + 3X^)E,{E, + uj) + {1 - X^){h ■ 1. + 1. ■ k)) , 
Lp'h;I = -UM^Ml ((1 + ?,\^)E,E, + (1 - ■ 1,) , (5) 

^ee,mixed ^ -UM^MI ((1 + 3\^)E,{2EI + E,UJ - h ■ k) + [l - \^)E,{2\, • 1, + 1, • k)) 

with nie the electron mass. In realizing the amplitudes from the Feynman rules we im- 
pose ej(fc) • A; = for each polarization state i of a real photon with momentum fc, noting 
^j=i2^i(^) ■ ^«(^) ~ To effect the subsequent photon polarization sums, however, 

we employ QED gauge invariance and make the replacement 2 — ^ ~9^^ 
throughout, without any supplemental conditions. 

Denoting the 0{a) correction to the amphtude by A^ioop the amended decay rate is 
determined by 

I-MT = l-Mtreel' + -Mtrec " -MLp + -^loop " A^Lc + 0{a^) . (6) 

The T-odd triple momenta correlation ^ = 1,^ ■ (le x k) in the decay rate can arise from 
the interference between the tree level amplitude A^tree and the anti-Hermitian parts of the 
one-loop corrections to it, so that ultimately the interference term ^spins(2Re(A^trccA^ioop)) 
contains terms linear in ^. Since we consider the decay and detection of unpolarized particles 
exclusively, X^spins l-^lx-odd indeed characterized by terms linear in ^. Evidently the 
induced asymmetry is suppressed by a factor of a = e^/47r ~ 1/137; explicit computation 
shows it to be much smaller still. 

Before turning to the computation of |Al|T_odd consider its relation to a measurable 
quantity. Following Ref. [H], we define a T-odd asymmetry A^, namely, 

where iV+ is defined as the total number of decay events with positive ^, and N_ is defined 
as the number of events with negative ^. Specifically, we compute 

where r± contains an integral of |A^p over the region of phase space with ^ < 0, respectively; 
the numerator is non-zero if and only if lA^lx-odd non-zero. Working to corrections of 
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0{e/M), the neutron radiative /3-decay rate F in the neutron rest frame is 



8M„ (27r) 



|ie|(i£;eo?ne'^(icc;(infc(int.^ j^j^^^ ^ — I - \m\ 



4M„, 



spins 



(9) 

where = M„ — le — lp — k and -E^ = M„ — Mp — Eg — u are fixed throughout. The precise 
form of r± depends on the concrete choice of coordinate system. Choosing the direction of 
the electron momentum Ig as the z direction and letting k and le fix the x-z plane, then under 
this specific choice ^ > corresponds to (pi, G [0,7r] and ^ < corresponds to cpu G [tt, 27r]. 
Thus we define 



r+(a;'"'") 



1 



16M2(27r)f 



ujduj I \\e\dEe I dxk I dXp 

J me J —c 



-1 
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spins 



(10) 



and 



16M2(27r) 



Lodoj 



%\dEe 



dxk / dxi 

c J-1 



2lT 



X 



spins 



Pp, Ey 



(11) 



where E^^^ 



Mn - Mp - CO, to' 



Mn — Mp — rrif., and u;™™ is determined by the 



threshold energy of the detector. In our computation of \Ai\ we set M„ = Mp = M in 
terms which would yield corrections beyond leading order in the recoil expansion. We limit 
the integration over Xk to the range [— c,c]; we discuss this as well as our choice for c in 
Sec.IVl 



III. COMPUTATION OF Espinsl-^lx odd IN LEADING ORDER 

To compute the T-odd pieces, we need to obtain the anti-Hermitian parts of the one- 
loop diagrams Im(A^ioop)- We do this by performing "Cutkosky cuts" [13], which means we 
simultaneously put intermediate particles in the loops on their mass shells in all physically 
allowed ways and then perform the relevant intermediate phase space integrals and spin 
sums. Graphically speaking, after imposing the cuts, the anti-Hermitian part of a one-loop 
diagram can be viewed as the product of two physical tree-level processes. We have 

Im(A^loop) = ^^jdpnY^ MfrM*n = ^ j ^PnYl ^fnMni , (12) 

Tl Sji Sfi 

where refers to the summation over all the possible cuts of the one-loop diagrams and 
/ dpn and Yls^ refer to the intermediate phase space integration and spin sums, respectively, 
for a cut which yields state n. The matrix elements M.ni and M. fn refer to the two tree- 
level diagrams after a physical cut. After excluding the physically unacceptable cuts, 14 
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FIG. 2: All two-particle cut contributions to n{pn) — >• piPp) + e~{le) + Veilp) + ^{k) which appear 
in 0{a) up to corrections of recoil order, using the syntax of Fig. [T] A "x" means that the 
intermediate particle has been put on its mass shell; two such symbols define the Cutkosky cut. 
The diagram enumeration is utilized in our calculation of the T-odd asymmetry. Note that the 
first number selects a particular Feynman diagram, and the second determines the particular two- 
particle physical cut in that diagram. 




FIG. 3: Compton scattering diagrams which appear in Im(7Wioop) for 7 — e cuts. We denote the 
two graphs by A^jjg(/g, A;', /g; ^) and A4'^^{1'^, k' ,le, k), respectively. The diagrams and amplitudes 
appropriate to 7 — p scattering follow from replacing electron with proton variables. 
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cut diagrams remain, and they are illustrated in Fig. [2j We evaluate them explicitly. The 
momenta labeled as k\ V^, and p'^ refer to momenta of intermediate particles. In performing 
the Cutkosky cuts, each particle in a pair of particles is put on its own mass shell. 

It is useful to categorize the cuts as per the sorts of processes involved. That is, M. /„ 
describes the manner in which select particles rescatter, so that we can have Compton scat- 
tering or electron-proton scattering, the latter with or without the emission of an additional 
photon. The family of diagrams given by (1), (2), (5.1), and (6.2) contain Compton scatter- 
ing from the electron, as illustrated in Fig. |3| whereas the family comprised of (3), (4), (7.2), 
and (8.3) contain Compton scattering from the proton. In these families A^/„ is captured 
by one of the following expressions: 

-M?e(C^',^e,A;): 
-M^e(^e,^',^e,A;) 

or 

M';piPp,k',pp,k) 

where e' = e{k'). Correspondingly, Ai 
amplitude, as per the form of A^oi and A^o2, with only some of the arguments changed. 
Technically we define a "family" to be those contributions to the T-odd correlation which 
cancel amongst themselves to yield zero when we replace e or e* by k or e' or e'* by k' as per 
the Ward-Takahashi identities. 

Furthermore, there is an intermediate phase space integral over the kinematically allowed 
phase space. For 7 — e scattering we have 

with Poe = le + k, whereas for 7 — p scattering we have 
with Pop = Pp + k. Collecting the pieces, we have 



Im(X,) = ^ j 




(19) 


Im(X.) = ^ j 


rfp^e ■^leO'e^ k' , k)Moi{l'^, k' , Pp) , 


(20) 


I.n(X,„) = ^^ 


j dp^e M^^{1'^, k', le, k)Mo2{l'e, k\pp) , 

^76 


(21) 




j dp^e J2 ■^7e(C' ^e, k)Mo2{l'e, k' , Pp) , 


(22) 




7 





= -e'uMe f''^lY^ /uM, (13) 

= eV(/e)^^^^Ae(0, (14) 

= -e'^UpiPp) '^^^ 2p -V ^ ^'^P^P'p'^ ' (-"-^^ 

= e^Up{Pp f^\p, . / ^ fMp'p) > (16) 
ni is given by the tree-level neutron radiative /3-decay 



FIG. 4: Diagrams which appear in Im(A^ioop) for e—p scattering with electron bremsstrahlung. 
We denote the two graphs by Mfp^{l'^,Pp,le,k,Pp) and Mfpy{l'^,Pp,le,k,Pp), respectively. The 
diagrams and amplitudes appropriate to proton bremsstrahlung follow from exchanging electron 
and proton variables. 



Pn 




FIG. 5: Contribution to n{pn) — )• piPp) + e {le) + Ve{lu) decay after Fig. [T] 



for the "7 — e" cuts, and 
Im(>l3) = 

Im(A^7.2) 
Im(A^8.3) 



^ j dp^p^M'^p{p'p,k',Pp,k)Mo2ile,k',p'p) 
^ j dp^P^M''^p{Pp,k',pp,k)Mo2ile,k',p'p) 
^ j dp^P^M'^p{Pp,k',pp,k)Moi{le,k',p'^ 
^ j dp^p^M''^pip'p,k',pp,k)Moi{le,k',p^ 



(23) 
(24) 
(25) 
(26) 



for the "7 — p" cuts. 

In addition to the families of Compton cuts, there are cuts in which A4 fn is determined 
by electron-proton scattering either with and without bremsstrahlung, and, correspondingly, 
Aini is determined by either nonradiative or radiative /3-decay. Referring to Fig. [2| we see 
for cuts in which the electron and proton scatter with bremsstrahlung that diagrams (5.2) 
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and (6.1) comprise the family associated with electron bremsstrahlung, as shown in Fig. |4| 
and (7.1) and (8.1) comprise the family associated with proton bremsstrahlung. In these 
families Ai jn is given by one of the following: 

-^:^,(/;,P;,/e,fc,Pp) = -e%e(/e) ^^^ ' ^f^ ^UeiQ ^p>^^l\^)2 ^M)l'' U,{p'^) , (27) 

M%,^{1:,p'^XXPv) = e%{pf^^^^^:^-^ , (29) 



or 



Mfp^{l'^,p'p,le,k,pp) = -e^np(pp)7^ ^^^^^^ _ ^ Up{Pp) ^^"^ ^^ UeiQl^UeiQ . 



(30) 



Moreover, M.ni is given by neutron /3-decay, as shown in Fig. |5} which, up to recoil-order 
corrections, reads: 

MMl'e.P',) = ^ne(/;)7p(l - l,)vMuMh'(^ - ^75)«n(Pn) • (31) 

Collecting the pieces, we have 

Im(A^5.2) = ^J dp,p,J2^epi(^e.P'pJe,k,Pp)MMl'e,P;), (32) 

Sep 

Im(A^6.l) = ^ j dp^p^^Mtp^{l'^,p'p,le,k,Pp)MY,K{l'e.P'p). (33) 



and 



Im(>l7.l) = ^/ ^Pep^E-^S7(^e,P;,^e,A;,p,)A^DK(/l,P;), (34) 

Sep 

Im(A^8.l) = ^ j dpep^^Mfp^{l'^,p'p,le,k,Pp)MBKil'e,P'p) (35) 



for the "e — p — 7" cuts. The last family of cuts is given by (6.3) and (8.2) in Fig. |2| In this 
case Ai fn is given by e — p scattering, and we have 

MepiCp'pJe^Pp) = -e'uMYnM ^'\^y np{pp)Yup{p'p) . (36) 

The corresponding J^ni is given by A^oi(C5^iPp) for (6.3) and A^o2(^ei ^iPp) (8-2)- We 
thus have: 

Im(A^6.3) = ^ j dpep^Mep{l'e,p'pJe,Pp)Moi{l'^,k,Pp) , (37) 

Sep 

Im(A^8.2) = ^ J dpepJ2^ep{l'e,P'pJe,Pp)Mo2{l'„k,p'p) (38) 
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for the "e — p" cuts. In these graphs the intermediate momenta satisfy l'^ + p'^ = le + Pp, 
so that the integral over the allowed phase space is slightly different from that in families 
with e — p scattering and bremsstrahlung. In particular, diagrams (6.3) and (8.2) are each 
infrared divergent when = 1^; this divergence cancels, however, as expected [22], once we 
construct A^^. 

The expressions we have collected complete the building blocks of the computation of 
the T-odd correlation in 0{a) up to recoil-order corrections. The spin-averaged T-odd 
correlation is 

y^T-odd = ^ E l-^lT-odd = I Yl('^Re{M,ree^lmMlJ) , (39) 

spins spins 

and we report the contributions to it family by family as each family represents a QED gauge 
invariant group of contributions. We employ a subscript system to identify the contributions 
in a straightforward way. Since the T-odd correlations are given by the interference of tree- 
level diagrams, which are numbered in Fig. [l] as (01) and (02), with one- loop level diagrams, 
which are numbered in Fig. |2] as (1), (2),..., and (8.3), we label, for example, the T-odd 
correlation from the tree diagram (01) and one-loop diagram (6.3) as |A^Pr|._Qj^ [6.3.01]. In 
computing the intermediate phase space integrals which enter these expresssions, we find 
that both vector and tensor structures appear in the intermediate momenta. We simplify 
such integrals using the Passarino-Veltman reduction [23] and present the details, as well 
as all needed integrals, in Appendix A. In Appendix B we report concrete expressions for 
the final gauge-invariant combinations of the various contributions to l-^Px-odd "which re- 
sult after performing the trace calculations and employing the formulae of Appendix A for 
the intermediate phase-space integrals. We work to leading order in the recoil expansion 
throughout. Judging by the structure of the resulting expressions one can see that some 
families, namely, the j—p family containing cuts (3) + (4) + (7.2) + (8.3), as well as the e—p—j 
family containing cuts (7.1)+(8.1), do not have leading-recoil-order contributions, whereas 
others do and need to be considered carefully. The computations necessary to determine 
Im(A^ioop) and the resulting T-odd interference term are involved, so that we employ the 
program "FORM" to compute analytic expressions for the traces [21]. We compute all of 
the diagrams with these methods as a check of our procedures - we verify that the expected 
cancellations do indeed occur. 

IV. RESULTS 

Before presenting our final results for the asymmetry, there are three important remarks 
to be made concerning our numerical evaluation of the integral of | Px-odd allowed 
phase space. First of all, we note that the contributions to the asymmetry from the e — p 
and e — p — 7 cuts dominate the final numerical result. The j — p contribution vanishes 
in leading order, whereas the 7 — e contribution partially cancels - the latter observation 
comes from our detailed numerical evaluation of the asymmetry. Second, we note that the 
contributions from the diagrams of the e — p cuts each contain an infrared divergence; we 
regulate this by inserting a fictitious photon mass m^. However, as we show in Appendix B, 
the infrared divergence cancels in the net contribution to the asymmetry from the e — p cuts. 
The remaining piece is thus finite and well-defined, and we can safely set to zero. Finally, 
we note it is most convenient to choose a restricted range in the 7 — e opening angle. As 
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TABLE I: T-odd asymmetry as a function of Wmin for neutron radiative /3-decay. 



^^min(,J-VleV j 






U.Ul 


i. ( D X iU 


-5 


0.05 


3.86 X IQ- 


-5 


0.1 


6.07 X IQ- 


-5 


0.2 


9.94 X IQ- 


-5 


0.3 


1.31 X IQ- 


-4 


0.4 


1.54 X IQ- 


-4 


0.5 


1.70 X IQ- 


-4 


0.6 


1.81 X IQ- 


-4 


0.7 


1.89 X IQ- 


-4 



one can see from the formulae in the Appendix A, the solutions to the Passarino-Veltman 
equations become invalid if the opening angle Oe-y between the outgoing electron and the 
photon is exactly equal to or to vr. There is no physical divergence. Rather, the spatial 
components of the vector and tensor equations to determine the relevant coefficients become 
degenerate at such a boundary. Potentially one could remove this difficulty by solving the 
equations for infinitesimal values of Oe-y or {6^^ — tt) and then interpolating the solutions to 
the needed 6^^ = and vr points. In our present work, we simply choose a restricted range 
Xk = cosOe^ € [—0.9,0.9], which spans the angular range over which the neutron radiative 
decay rate is largest [23] . 

We can now present our results for A^. Noting Eq. ([8|, we see that Eqs. and ([5| 
share a common factor of e'^g'yG'jpM'^ /2, making independent of the decaying particle's 
mass in leading order in the recoil expansion. As can be seen explicitly in Appendix B, 
all of the contributions to lA^Px-odd found to be proportional to (1 — A^), so that 
the resulting asymmetry goes as (1 — A^)/(l + 3A^), up to small corrections, in this limit. 
The universal dependence on A in lA^Px-odd stems from the special nature of the T-odd 
correlation. It is a real triple product in momenta arising from the interference of a tree-level 
diagram with an imaginary part of an one-loop diagram after summing over the particles' 
spins. To leading order in M, the only surviving contribution is obtained from the product 
of the symmetric part of the lepton tensor, which is determined by a trace containing 75, 
namely, l^e"'^^^ + l^e"^'^^ — g^'^ly ^e^^'^^ , where a, /3, and 7 refer to photon or lepton indices, 
with the symmetric part of the hadron tensor. The latter is proportional to (1 -|- \^)ppPs — 
X^M'^QpS, where p is a baryon momentum and = M^. As one can easily check, this special 
combination generates an overall (1 — A^) coefficient; the remaining (1 + A^) term cannot 
be of leading order once the photon spin sum is effected. We use A = —1.2701 ± 0.0025 [26] 
in our numerical evaluation. For definiteness, the remaining input parameters we employ 
are rrie = 0.510999 MeV, M„ = 939.565 MeV, Mp = 938.272 MeV, and a'^ = 137.0360 - 
these quantities can be regarded as exact for our current purpose [SE]- We show our results 
for the T-odd asymmetry in neutron radiative /3-decay in Table [T] and Fig. |6} We see that 
the asymmetry is rather smaller than a. We recall that the radiative /9-decay rate grows 
as logWrnin as Wmin — )■ 0, whcreas the l-M^T-odd tends to zero in that limit. Consequently 
the small values of the asymmetry as Wmin — )■ is reflective by the growth in the decay rate 
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»min(MeV) 

FIG. 6: The asymmetry versus the smallest detectable photon energy Wmin in neutron radiative 
/3-decay. 



itself. 

Generally is determined by an interplay between A and the energetics of the decay, 
along with the value of Wmin. The universal (1 — A^)/(l + 3A^) behavior in A makes for a 
rich pattern. If, for some nucleus, the associated value of A were significantly different from 
unity, the T-odd effect could be considerably amplified, whereas for if A ~ 1, the T-odd 
effect could be substantially reduced, facilitating from this perspective at least the search for 
physics beyond the SM. Interestingly a "quenching" of the Gamow- Teller strength in nuclei 
in relation to shell-model predictions is experimentally established [271 EH] ^ it derives from 
the presence of many-body correlations in the nucleus ^29j. As an illustration, we consider 
the process ^^Ne — )■ ^^F + + Ue + j. We map our formalism to this new case by using the 
parameters M^nc = 17700.13744 MeV and Mwp = 17696.898613 MeV [30], as well as gf = 
0.928, as determined by Refs. [311 E2] with Ref. [S3] for a translation from the conventions 
of those references to g'^. The axial- vector coupling can be determined in /3-decay from 
either the measured decay rate [211 E2] or the /3-asymmetry [51] , and the determined values 
are in agreement, though the former method yields a more precise result. To determine 
A^^ = gf/gv, we note gv = {l + A^y/^V^A, where we use Kd = 0.97425 ± 0.00022 [26] and 
= 0.0240 [35], though further refinements of the latter are known in principl e \36\ l37] . 
We employ A*^^ = 0.941 and present the asymmetries as a function of a;™™ in Table jn} As one 
can see from the numerical results for fixed Wmm, the asymmetry in ^^Ne radiative /3-decay 
is in fact significantly smaller than that in the neutron case. 

In our paper, we compute the 0{a) contribution to the T-odd asymmetry, keeping only 
the leading terms in the recoil expansion. The accuracy of our calculation is limited by 
the uncertainties in the input parameters we employ, as well as by the numerical size of the 
neglected recoil-order contributions. Crudely we expect the latter to be reduced with respect 
to the leading-order contribution by a factor of 0{E^^^/M) ~ 1 ■ 10^'^. The uncertainty 
in A^^ may well be the dominant error in some nuclei. Nevertheless, we can conveniently 
check the rough size of the recoil-order contributions in the neutron case by replacing the 
vertex ^^{gv — gAlb) iii the tree- level amplitude with the weak magnetism contribution, 
—ia^^qyF2{q^)/{2M), where F2(fi)/gv = i^v = 3.706, the isovector magnetic moment of 
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TABLE II: T-odd asymmetry as a function of uj^i^ for Ne radiative /3-decay. 



Wmin(MeV) 








0.01 


-2.39 X 


10" 


-6 


0.05 


-3.96 X 


10- 


-6 


0.1 


-5.27 X 


10- 


-6 


0.2 


-7.31 X 


10- 


-6 


0.3 


-8.97 X 


10- 


-6 


0.4 


-1.04 X 


10- 


-5 


0.5 


-1.16 X 


10- 


-5 


0.6 


-1.26 X 


10- 


-5 


0.7 


-1.34 X 


10- 


-5 



the nucleoli. The interference of the resulting recoil-order contribution with the tree-level 
amplitude yields upon explicit calculation a contribution to which is no larger than 
~ 4 ■ 10-^ for Wmin = 0.3 MeV. 

V. SUMMARY 

In this paper, we have computed the T-odd correlation in neutron radiative /3-decay aris- 
ing from SM physics. The T-odd correlation is characterised by the kinematical variable 
^ = 1,^ ■ (le X k); consequently, it is spin-independent - and thus fundamentally different 
from a permanent EDM. The mimicking T-odd correlation arises from the presence of elec- 
tromagnetic final-state interactions when the intermediate particles are each put on its own 
mass shell. Taking the neutron prototype, we have computed the leading-order 

result, which is of 0{a, (e/M)°), to the T-odd asymmetry exactly. In particular, our de- 
tailed analysis shows that the resulting T-odd asymmetry is controlled by (1 — A^), so that 
vanishes as A — )■ 1, suggesting that radiative /3-decay studies in other systems could be 
employed to good effect. Indeed our procedure is completely general, so that the resulting 
formulae can be readily translated to hyperon and nuclear radiative /3-decays as well. We 
have explicitly considered the example of ^^Ne radiative /3-decay as a concrete illustration of 
the possibilities. Its lifetime is much shorter than the neutron's, making experiments more 
practical, and it should be possible to study such decays in a trapped atom experiment 
It also possesses a value of A''^ rather close to unity. The combined effect of these obser- 
vations is a T-odd asymmetry roughly an order of magnitude smaller than that in neutron 
/3-decay. In this case A^^ presents a much reduced background against which to search for 
new sources of CP-violation. 
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Appendix A: Intermediate Phase Space Integrals 

The computation of the imaginary parts of the loop diagrams requires an integration 
over the allowed phase space of the intermediate momenta as fixed by the momenta of the 
final-state particles and energy-momentum conservation. In this Appendix we report the 
integrals which appear in the diagrams of Fig. |2] and label them as per the diagrams in that 
figure. For diagrams with cuts which yield Compton scattering from electrons our results 
can be compared to, and agree with, those of Refs. [H] and [12] • In what follows we report 
the integrals which arise from 7 — e cuts: (1), (2), (5.1), and (6.2), and then the integrals 
which arise from the cutting of electron and proton lines to generate physical ep — )■ epy 
scattering, namely, (5.2) and (6.1), and ep — )■ ep scattering, (6.3) and (8.2). The integrals 
associated with the rest of the cuts in Fig. |2] are not given explicitly because they do not 
contribute in leading order in the recoil expansion, as we note in the main body of the text. 
Nevertheless, we note the relationships between these integrals which appear in the large 
Mp limit in order to make the cancellations associated with these terms transparent. 

From diagram (1), defining Pqe = h + k, we have 



'76 

,2 



as well as 



with 



= / dp^.k'f^ = a,P^^ (A2) 



2 



From diagram (2) we have 



TT / ml 



1 TT . 



We apply the Passarino-Veltman reduction method to compute integrals which contain ad- 
ditional powers of the intermediate momenta |i23|. That is, writing 



= I dp^ej^ = a,li: + h,Pi;, , (A4) 



the values of 02 and 62 are fixed by the solution of the set of equations 

= a2ml + 62/e ■ -Poe , 

If. ■ kJ2 = ^2^6 ■ -foe + ^2-Poe • 
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Moreover, 

Lf = j dpj^^ = c,g^^ + d,l^Jl + e^Po'en^; + U^P^ + P^Z) , (A5) 
where C2, (i2, 62, and /2 are given by the solution of the set of equations 

= 4c2 + d2ml + e^Pl + 2/2/e • Poe , 

= C2 + (i2mg + /2/e • -Poe , 
ai = e2le ■ Poe + /2"^e > 
/e ■ /c62 = C2 + e2Poe + /a^e " Poe ■ 

For integrals which depend on Mp wc report their form in the large Mp limit for subsequent 
use. Note that M rather than Mp appears in the limiting form because the n — p mass 
difference itself is of higher order in the recoil expansion. Prom diagram (5.1) we have 

7 1 ^1 fPp-Poe + Ioe\ 



with loe = J{Pp- PoeY - M^P^, noting 



7 ^ /^e + W+ |k + le| , 



2M k + le \Ee + uj -\k + U 



as Mp ^ 00. In addition 



r u'n 

Ki = J dPye^ = a,.,p; + 65.1P0': , (A8) 
where as.i and 65.1 are given by the solution of the set of equations 

Jl = a5,iMp + h.lPp ■ Poe , 
le ■ kJ^^i = a^.lPp • Poe + ^5.lPoe ■ 

In the large Mp limit 65.1 ~ 1/M and as.i ~ 1/M^. We postpone discussion of the integrals 
from diagrams (5.2) and (6.1) to consider the integrals from the remaining diagrams with 
Compton cuts. From diagram (6.2) we have 



with /e = \ Kpp ■ leY — M^ml and 



~ , log ( ) . (AlO) 



2M\L\k-L °\E,-\l 



as Mp — > 00. In addition 
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where 06.2, &6.2, and ce.2 are given by the solution to the set of equations 

J2 = aQ,2Pp ■ Poe + h.2Pp ■ le + CqJiM^ , 
J5.I = 06.2^6 • Poe + h.^ml + C6.2Pp • h , 
le ■ kJe.2 — Q6.2-Poe + ^6.2^e • Poe + CQ.2Pp ' Poe ■ 

and in the large Mp limit aQ,2, &6.2 ~ 1/M and cq,2 ~ 1/M^. Finally 

W2-J ^^7e(;^.^,)(^^.^,) 

= (1^.29'" + e6.2pX + + 96.2Po'ePoe + /i6.2(pj^/r + ^e^Pp) 

+^e.2{p'pPoe + Op) + fc6.2(/e^Po'e + PoTe) , (Al2) 

where the coefficients which appear are given by the solution to set of the equations 

44.2 + 66.2^^ + f6.2ml + fl'6.2-Poe + 2/i6.2Pp ' k + '^i&.2Pp ' Poe + 2/^6.2^6 " -Poe = , 

de.2 + h.2ml + ^6.2Pp • U + ^6.2^e • -Poe = , 

e6.2Pp • le + ^6.2"^e + ^6.2^6 ' -Poe = 05.1 , 

5'6.2^e • -Poe + «6.2Pp ' U + A;6.2"T.e = 65.I , 

4.2 + 66.2-^^ + /l6.2Pp • + «6.2Pp " pQe = , 

5'6.2Pp • Pbe + i&.2Ml + A;6.2Pp -le^h, 

dQ.2 + 5'6.2p0e "I" '^6.2Pp ' Poe + ^6.2^e " pQe — le ' ^06.2 ■ 

Note that the equations have been chosen to yield a self-consistent solution for the six 
coefficients. 

The integrals associated with the 7 — p cuts can be found if necessary by replacing the 
intermediate momentum by p'p as well as /g by Pp in the 7 — e integrals we have provided. 
Specifically we note 

where Pop = Pp + k, and 



so that 



as Mp — >■ 00. Moreover, 



and 



whereas 



J4 = / dp^p ^ (A13) 
J Pp-k' 

J7.2 = J (^P-ypf—^ 



7.2 = J (^PyPfry = «7.2/e^ + b7.2p; , (A16) 
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and 



so that 



68 3 ^a7 2 + C'f^) ; C8 3 —b72 + o(^] (A19) 



as Mp — oo. 



The integrals in the remaining diagrams of Fig. |2] arise from cutting the electron and 
proton lines to generate physical ep — )■ ep7 or ep — )■ ep scattering. The intermediate phase 
space integrals in these cases are more complicated than those associated with the Compton 
cuts; fortunately, closed-form expressions for the integrals in the large Mp limit suffice to 
leading order in the recoil expansion. With Pg = Pp + + k, we note for diagram (5.2) 



(P2 _ M2 + m2)2 - AP^ml ~ —^/iE, + uy-ml (A20) 



2p2 V -p ■ M 



as Mp — )■ oo. Moreover, 

Jb.2 = I dpepy—^ — (A21) 



and 

TT 



J, , - log I ml + h-k-iE^ + uY + ^iE^ + u)^-ml\l^ 

4M|le + k| ^\ml + i^.k- {E, + uy - ^{Ee + u)^ - ml % + k| ' 



as Mp — )► oo. In addition. 



^5.2 = / dp^p^—^— = a,.2Po^, + c,.2P^p , (A23) 
J \Pp Pp) 



where 05.2 and C5.2 are given by the solution to 

2 



ml + le- /c) J5.2 - "TT = «5.2-Poe + C5.2Pp " Poe , 



1 

Pp ■ PoeA2 + = a'5.2Pp ' Poe + £5.2^^ 



p ) 



SO that in the large Mp limit 05.2 ~ 1/M and C5.2 ~ Turning to the integrals from 

diagram (6.1) we have 



so that as Mp — )■ 00 



17 



as well as 



where as Mp oo 



^''■'^ ^^^^^ 

~ 2{ml + h ■ k)hA - 2/5.2 - 2/6.1 (A27) 



with 



^6.i = ^rr:^ I ViE, + uy-ml + — j log I — j===== 

(A28) 



Moreover, 



so that as Mp — >■ 00 



where 

A± = + • A; - (^e + t^)' ± |le + k| V(^e + - ml (A31) 

and 

B± = |le|'(/e ■ kf - {u^ml - E,u{le " k)) A± 

+ |le|(/e • k) ({E, + a;)a;|le + k| ^ (^^' + le • k) ^(Ee + ujf - mi) . (A32) 

In addition, 

Ki = / dpepjjjr-rxf = "6.1^^ + ^'e.iA;'^ + ce.ip^ , (A33) 

^ l^e ■ '^J IPp ~ Pp) 

where the undetermined coefficients are fixed by the solution to 

J5.2 = ae.ile ■ k + Ce,iPp ■ k , 
{ml + le ■ k)jQ,i — = aQ,i{ml + /g ■ k) + h^Je • k + Ce.iPp • Poe , 

so that in the large Mp limit Oe.i, fce.i ~ 1/-^ a-nd Ce.i ~ 1/M^. Also 
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where the undetermined coefficients are fixed by the solution to 

44.1 + ee.iMp^ + fe.iml + 2hQ,iPp ■ k + 2i6.iPp ■ k + 2A;6.i/e ■ k = Je.i , 
4.1 + ee.iM^ + h.iPp ■ le + ie.iPp ■ k = Pp ■ PoeCe.i , 



96.iPp ■ k + ie.iMp + fce.iPp ■ h = Pp ■ -Poe&e.i 
r2 



ee.iPp • + hG,ile ■ k = C5.2 , 
fe.ile ■ k + hQ,iPp ■ k = 05.2 , 

de.lPoe + e-e.iiPp ■ Poef + fe.lih ■ Poe)^ + ge.lih ■ kf + 2/l6.lPp • Poek ■ Poe 



+2i6.iPp • Poek ■ k + 2kQ,ile ■ Poeh ' k = {ml + le ■ kfJeA - {ml + le ■ k)h,i + -|i . 
For the remaining e — p — 7 cuts we have 

K^.i = j dPe^Jf±J^ = + br.ip^ , (A35) 

''■'^l'''^{p'p-k)li'e-le? ^^''^ 
^8^1 = j dpep ^ _^ )2 = «8-l^e + ^S-l^'' + Cg.lP; , (A37) 



and 



and 



whereas 



so that 
as Mp — >■ 00. 

The integrals for the e — p cuts follow from those we have just analyzed under the 
replacement of Pq with Pg = Zg + Pp- In this case, however, there is an added complication 
because the integrals become infrared divergent when p'^ = Pp. This divergence cancels once 
we construct an observable quantity; nevertheless, we regulate the integrals as they stand 
by adding a fictitious photon mass - this will allow us to track the infrared divergences 
through the course of the calculation, so that we can demonstrate the divergence cancellation 
manifestly. In what follows we set to zero in all terms which are finite in the 
limit. We have 
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as Mp — > oo. In addition, 



log 



7 



as Mp — )■ oo. Thus we see that Js.2 vanishes in this hmit save for the infrared divergent 
piece, which we define as Jg'J. In addition, 

r 1 I'l^ 

= / <ll>..^,J^^p^^^ = asJ^ + + ^s.!^ . (A41) 

The coefficients are given by the solution to 

mlJ8.2 - = a8.2ml + bsJe ■ k + C8.2Pp ■ le , 

Pp ■ leJs.2 + ^^8.2 = as.2Pp ■ le + h.2Pp ' k + 08.2-^"^ , 



where 



2 

{le + Pp) ■ kJs.2 - -^8.2 = «8.2^e • k + C8.2Pp • k , 



1 



p'p-k ' y ^''{i^p-Pp) 



8.2 = I dpep—r ; /8.2 = / dpepj-, ^ ■ (A42) 



2 — 



In the large Mp limit we note that 



K^.2 - ^^2 (A43) 



so that 68.2 ~ 0, and we need only solve 

"^e^8.2 - = a8.2ml + C8.2ME^ , 

EeJ8.2 = 08.2^6 + C8.2M (A44) 

to determine the leading-order expressions for 08.2 and C8.2- We can track the infrared 
divergence in J8.2 in 08.2 and C8.2 by solving these equations with 78.2 = and J8.2 = ^s.'J; 
which yields af^ ~ J|'J and Cg'^^ ~ in leading order. 
The integrals from diagram (6.3) are 



'6.3 



log 



2uM V-Ee-IIel 



as Mp — >■ 00 and 



- j dP^v\^ (A46) 
~ 2mlh,2, - 2/6.3 
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with 

? TT ( _ , ^ , , /-Ee + 11 



'6.3 



(^(£;,2-£;e|le|c0S^e)l0g(|^^) +211^^008^^) (A47) 



IMuj 

as Mp ^ oo. We define k • Ig = |k||le| cos6'e. Moreover, 

~ loe — h loe — - — 

^\\\{U-k)M \ ^4|1,|2^ ^{U-hf 

as Mp — )■ oo. In this case we see that Je.s has both infrared finite and divergent pieces in 
the Mp — >■ oo hmit - the latter we define as Jg J. Finally 

where the undetermined coefficients are fixed by the solution to 



^6.3 - = Oe.S'^e + ^6.3^e " + C6.3Pp " le , 



Also 



/■ 1 I'^l"" 

= J ^^-^ (Z;.fe)fa-;;)2-^2 = ^6.3^'^'^ + e6.3Pp%^ + hJX + 9e.3k^k^ 

+h,.^Wl + I'^vl) + ^6.3(P^^^'^ + + kUKk"" + A;'^/:) , 

where the undetermined coefficients are fixed by the solution to 

4(^6.3 + e6.3-^p + h.zml + 2/i6.3Pp • U + '^k.zPp ■ k + 2kQ.^le ■ k = Je.a , 
d&.z + ee.sMp + he.sPp • U + i6.3Pp ■ k ^ Pp • IgCe.s , 
g&.zPp ■ k + ie.sM^ + ke_3Pp ■ h ^ Pp ■ kKz , 

h.zPp ■ le + ^6.3-/W"p + /Ce.sPp -k^pp- leaQ.2. , 
e6.3Pp • ^ + ^6.3^e - k^Pp- kCs.2 , 

faJe ■ k + he.3Pp ■ k ^ Pp ■ ka8.2 , 

4.3"^e + 66.3 (Pp -kT + f 6.3^1 + ge.sik • k^ + 2/i6.3Pp " + 2iG,3Pp ■ IJe ' k 

J' 

+2kG_sm% • k = mlJe,3 - mg/6.3 + ■ 

We can track the infrared divergence in Jg.s in the solutions for the vector and tensor 
coefficients by solving the equations in the large Mp limit with /q.s ~ /g 3 ~ and Je.s ~ Jq^s, 
with as.2 ~ Os.'J) which yields Og^J ~ /g'3 ~ JgJ with all other coefficients zero in this limit. 
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Appendix B: 1-^4 Prp^^^ in Leading Order 

In what follows we report the contributions to the T-odd correlation in 0{a) up to 
corrections of recoil order. We organize the results as per the various gauge-invariant families 
we describe in the main body of the text, employing the subscript convention which follows 
the labeling in Figs. [T] and |2} We use the integrals and Passarino-Veltman coefficients defined 
in Appendix A. The result for the 7 — e family is 



l-^IV-odd [1-01 + 1-02 + 2.01 + 2.02 + 5.1.01 + 5.1.02 + 6.2.01 + 6.2.02] 

mj ^ ml ml ^ MEe, ME^ ME^, 2MEe 

+ -J ; h + z ; a2 - J2 H ^6.2 H 96.2 06.2 ^6.2 

le-KUJ le-KUJ le-KUJ UJ UJ UJ LO 

ME,, MEe ME, MEe ME, ME, M^ . 
-J6.2 + -J ; — O5.1 ~ ; — ^5.1 ~ 777 r5'6.2 + r/6.2 + r'^6.2 H "^6.2 



CO ' le ■ kco ' le ■ kcu ' 21, ■ k ' 21, ■ k ' I, ■ k ' uj 
M2 M'^E, M^E,. M2 



2uj ^'^ 1, ■ k ^'"^ ~^ I, ■ k ^'"^^ 2l,-ku ^'^ ~^ 21, ■ k ^'^ 
The result for the 7 — p family is 



l-^IV-odd [3-01 + 3.02 + 4.01 + 4.02 + 7.2.01 + 7.2.02 + 8.3.01 + 8.3.02] 

9 9 ^9 .9.f E, E, 1 ME,, M 

\le-kuj I, - kio I, ■ koo^ I, ■ k 00 

^ME,^ ME, J ^ ^ ^ ^ J c ■ 

= + C(M) , 



where we employ Eqs. (A14) and (A19) to determine that the contribution to this family 
vanishes in leading order in M. The results for the e — p — 'j families are 



l-^pT-odd [5-2.01 + 5.2.02 + 6.1.01 + 6.1.02] 

99^9 ,oY2Ee, 2M M 2E, 1 2ml, 

= -a^glGliUM^l - ^fce.i + ze.i ce.i - ^-05.2 - ^K.i 

\ u u u 1, ■ kuj 1, ■ k 

ml ml Ml 2ME, 2ME, 2ME, M^ 

if rv if rv if rv UJ if n, if rv if rv if rv 



and 



l-^IV-odd [''•l-Ol + 7.1.02 + 8.1.01 + 8.1.02] 

2E, 

-J ; 0-7.1 + -, -, "7.1 r ^"7.1 "7.1 

Ip ■ kuj I, ■ kuj Ip ■ k 00 



9 9 ^9 ,9x/ 2E, M , 2ME, 2M, 

-a^g^GUQAM^l - A^) | ^^a7.i + ^^&7.i - -^«7.i - — &7.1 



= + 0{M) , 
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where we employ Eq. (A38) to determine that the contribution to this family vanishes in 
leading order in M. We emphasize that the contributions which vanish do so simply to the 
order of the recoil expansion in which we work. Finally, the result for the e — p family is 



T-odd 



[6.3.01 + 6.3.02 + 8.2.01 + 8.2.02] 



2ml 



2E, 



2E, 



2M 



M 



Ip ' k 



^6.3 



2ml 



06.3 



m; 
Ip ' k 



Je.3 + 



2MEe 
Ip ' k 



;.2 + 



u 

2ME, 
Ip ' k 



-«6.3 



2ME^ 
Ip ' k 



-«6.3 



-H.3 



-C6.3 



Ip ' k 



C8.2 



/g ■ k 



If, ■ k 



66.3 



6.3 



From Appendix A we note that a^-^ ~ .y--„- ~ 



J^'^ with 



i'2 ~ 4^ ~ ih ■ k)Ji^/{Muj) and afl ^ - ,6.3 
all other coefficients zero in leading order in the recoil expansion. Thus we see explicitly 
that the infrared divergence really does cancel in (9(M^). 
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